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ABSTRACT 

In  this  paper  we  explore  a  multiresolution  approach 
to  solving  one  dimensional  inverse  problems.  The  ap¬ 
proach  we  take  is  motivated  by  the  work  of  Chou, 
Golden,  and  Willsky  [1]  and  Beylkin,  Coifman,  and 
Rokhlin  [2].  Specifically,  we  consider  inverse  prob¬ 
lems  described  by  that  class  of  operators  which  are 
made  sparse  under  the  action  of  the  wavelet  trans¬ 
form  [2].  Moreover,  statistically- based  inversion  pro¬ 
cedures  utilizing  multiscale  a  priori  stochastic  models 
are  considered  [1].  As  a  concrete  example,  we  examine 
a  deconvolution  problem  arising  in  wellbore  induction 
measurement  of  conductivity. 


1  INTRODUCTION 

In  this  work  we  consider  the  utility  of  applying  mul¬ 
tiresolution  methods  to  the  solution  of  linear  inverse 
problems.  We  take  this  approach  above  all  others 
specifically  because  many  of  the  issues  arising  in  the 
solution  of  linear  inverse  problems  are  exactly  those 
which  are  handled  well  by  multiresolution  methods. 
First,  there  exist  inverse  problems  for  which  the  phe¬ 
nomena  under  consideration  may  be  fractal  in  nature 
or  characterized  in  an  interesting  way  by  their  behav¬ 
ior  on  a  multitude  of  temporal  or  spatial  scales.  Alter¬ 
natively,  the  data  upon  which  the  inversion  is  based 
may  exhibit  muitiscale  characteristics  as  is  the  case 
in  many  sensor  fusion  problems  [1].  Finally,  even  if 
neither  the  phenomena  nor  the  data  is  inherently  mul¬ 
tiscale,  recent  work  by  Beylkin,  Coifman,  and  Rokhlin 
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[2]  suggests  that  multiscale  methods  are  a  useful  tool 
for  purposes  of  computational  efficiency. 

We  consider  inverse  problems  described  by  that  class 
of  operators  which  are  made  sparse  under  the  action 
of  the  wavelet  transform  [2].  Thus,  the  inversion  pro¬ 
cedure  is  carried  out  in  the  wavelet  transform  domain. 
As  with  many  inverse  problems,  difficulties  associated 
with  ill-posedness  and  ill-conditioning  are  overcome 
via  regularization.  Because  the  inversion  is  performed 
in  scale-space,  we  employ  a  statistically-based  regular¬ 
ization  technique  utilizing  multiscale  a  priori  models 
similar  to  those  explored  by  Chou  and  Willsky  in  [1]. 

The  vehicle  for  our  work  is  an  inverse  problem  aris¬ 
ing  in  the  area  of  geological  exploration.  Specifically, 
we  are  interested  in  determining  the  conductivity  pro¬ 
file  about  a  wellbore  based  upon  a  suite  of  induction 
measurements  each  of  which  contains  conductivity  in¬ 
formation  at  different  spatial  resolutions.  This  prob¬ 
lem  is  of  particular  interest  because  it  fulfills  all  three 
criteria  for  the  application  of  multiresolution  tech¬ 
niques  to  inverse  problems.  That  is,  conductivity  vari¬ 
ations  are  typically  fractal  in  nature,  the  data  conveys 
information  on  a  variety  of  scales,  and  the  linear  oper¬ 
ators  describing  the  measurements  process  are  nearly 
diagonalized  by  the  wavelet  transform. 

2  INDUCTION  LOGGING 

Induction  logging  is  a  common  technique  used  in  ge¬ 
ological  exploration  for  determining  the  conductivity 
profile  of  a  rock  formation  surrounding  a  wellbore.  As 
shown  in  Figure  1,  an  induction  log  is  recorded  as  a 
transmitter  and  receiver  coil  pair  are  pulled  through 
the  formation.  A  current  in  the  transmitter  induces 
sympathetic  currents  in  the  rock.  In  turn,  these  for¬ 
mation  currents  create  a  secondary  magnetic  field  re¬ 
sulting  in  measurable  current  at  the  receiver  [3].  The 
problem  is  to  determine  the  conductivity  of  the  for¬ 
mation  given  the  current  signal  at  the  receiver.  More 
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sophisticated  induction  logging  techniques  involve  the 
use  of  multiple  transmitter  and  receiver  pairs  all  mea¬ 
suring  the  same  rock  structure  with  different  pairs  sen¬ 
sitive  to  different  areas  about  the  borehole.  As  the 
length  of  each  log  becomes  long  and  the  number  of 
antenna  pairs  grows,  there  exists  a  problem  in  terms 
of  performing  the  inversion  in  reasonable  time  and  at 
reasonable  cost  while  still  making  good  use  of  the  vast 
quantity  of  available  data. 

In  general,  the  conductivity  in  the  formation  is  re¬ 
lated  in  a  nonlinear  way  to  the  currents  observed  in  the 
formation.  Commonly,  the  problem  is  approached  by 
linearizing  this  relationship  about  some  a  priori  back¬ 
ground  conductivity  profile  [4].  Moreover,  it  is  often 
assumed  that  the  formation  is  azimuthally  symmet¬ 
ric  about  the  borehole.  In  this  paper,  we  reduce  the 
problem  to  the  one  dimensional  case  by  making  the 
further  assumption  that  the  conductivity  only  varies 
vertically.  In  this  case,  it  can  be  shown  [4]  that  the 
observations  at  the  ith  receiver,  gi(z)  are  related  to 
the  conductivity  field,  /(z),  through  the  convolution 
equation 

9i(z)  =  j  Ti(z  -  z')f(z')dz '  =  Ttf  (1) 

Where  Ti(z)  is  a  linearized  kernel  derived  from 
Maxwell’s  equations  and  dependent  upon  the  back¬ 
ground  conductivity  profile.  Two  examples  of  the 
Ti(z)  from  Schlumberger’s  Array  Induction  ToolrAf 
are  given  in  Figures  2  demonstrating  how  information 
at  different  resolutions  may  be  available  to  the  inver¬ 
sion  routine. 


Figure  1:  Geometry  of  Induction  Logging 


3  MULTIRESOLUTION  PROBLEM 
STATEMENT 

In  [2],  Beylkin,  Coifman  and  Rokhlin  proved  that 
a  large  class  of  operators  are  made  sparse  under  the 
action  of  the  wavelet  transform  and  derived  a  fast  algo¬ 
rithm  for  evaluating  integrals  containing  these  kernels. 
The  kernels  we  consider  in  conjunction  with  the  induc¬ 
tion  logging  problem  are  also  of  the  class  discussed  in 
\?}.  Thus,  we  are  able  to  exploit  the  sparse  and  nearly 
diagonal  structure  of  these  operators  in  the  wavelet 
domain  to  construct  their  inverses  using  computation¬ 
ally  efficient  algorithms. 

In  this  work,  we  adopt  the  following  notation.  In  one 
dimension,  the  scaling  function  coefficient  for  a  func¬ 
tion  f(z)  corresponding  to  the  basis  function  at  scale 
m  and  shift  k ,  is  given  by  f(m;k)  while  the 

wavelet  coefficient  associated  with  the  wavelet 
is  <f>(m\k).  If  the  wavelet  transform  of  a  function  is 
to  be  considered  only  over  a  finite  number  of  scales 


Figure  2:  Two  Induction  Kernels 
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them  the  coarsest  scale  is  denoted  L  and  the  finest  by 
M .  Additionally,  in  many  cases,  it  is  useful  to  con¬ 
sider  the  collection  of  scaling  or  wavelet  coefficients 
at  a  particular  scale  as  a  single  entity.  Hence  we 
adopt  the  notation  <j>(m)  =  {<f>(m;k)\k  £  Z}  and 
/(m)  =  {/(m;le)|A  £  Z}.  When  considering  the 
wavelet  transform  of  a  function  over  a  finite  number 
of  scales:  L  <  m  <  M,  we  use  the  following  construct: 

W(/)  =4>=  [<t>T{M  -!)•••  <1>T{L)  fT(L)f  (2) 

where  W(/)  denotes  the  operation  of  taking  the 
wavelet  transform  of  the  function  /.  In  two  dimensions 
we  consider  only  tensor  product  wavelets  [2].  Here, 
the  fundamental  grouping  of  wavelet  coefficients  for 
the  transform  of  T(z,z')  is  the  standard  form  repre¬ 
sentation,  WT<W*  =  ©,  as  discussed  in  [2]. 

Using  our  notation,  Equation  1  may  be  rewritten  as 

T  =  W(!7)  =  [WTiW][W(f)\  =  e«t>.  (3) 

The  object  of  the  inversion  is  to  construct  and  apply 
0  r 1  so  that  we  may  obtain  <f>  given  the  wavelet  trans¬ 
form  of  the  function  g(z).  Then,  f(z)  is  recovered  via 
the  inverse  wavelet  transform. 

An  example  of  the  structure  of  ©  for  the  conductiv¬ 
ity  deconvolution  problem  is  given  in  Figure  3  for  the 
convolution  kernel  corresponding  to  Z\(z)  and  using 
the  Haar  wavelet  for  the  transform.  Following  [5],  we 
expect  the  inverse  of  ©x  to  be  sparse  and  constructible 
in  a  computationally  efficient  manner.  Unlike,  [5]  how¬ 
ever,  we  propose  a  stochastic  inversion  scheme  based 
upon  a  multiscale  a  priori  probabilistic  model  for  the 
wavelet  transform  of  f(z). 

4  A  STOCHASTIC  REGULARIZATION 
APPROACH 

Two  common  difficulties  associated  with  inversion 
problems  of  the  form  considered  in  Equation  1  are  ill- 
posedness  and  ill-conditioning  of  the  discretize  system. 
The  solution  to  these  issues  usually  comes  in  the  form 
of  a  regularization  procedure  designed  to  constrain  the 
original  problem  so  as  to  make  it  well-posed  and  to  re¬ 
move  the  numerical  difficulties  associated  with  poorly 
conditioned  linear  systems.  In  this  work,  we  choose  a 
probabilistically  based  regularization  approach. 

4.1  The  Regularization 

The  regularization  is  accomplished  in  two  steps. 
First,  we  assume  that  the  measurements  are  corrupted 
by  an  additive,  white  Gaussian  noise  process,  n(z),  so 
that  the  input  to  the  inversion  algorithm,  t/(z),  is  given 
by  _  ' 

y«(*)  =  y«(z)  +  nj(z)  (4) 

=>  W(yj)  =  rn  =  &i(f>  +  V.  (5) 
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Figure  3:  Grayscale  Image  of  |©x|  Using  Haar  Wavelet 

where  v  =  W(n).  Second,  we  construct  a  prior  model 
for  the  function  f(z).  Because  we  are  interested  in 
performing  the  inversion  in  the  wavelet  transform  do¬ 
main,  a  natural  choice  for  a  prior  is  a  multiscale  model 
of  the  form  considered  [1].  Here  we  take  the  model  as 

/(m  +  1)  =  £f*(m)/(m)  +  (6) 

<t>(m)  ~  A/- (0,  R(m))  (7) 

E[<t>{m)4>(n)*}  =  R(m)6(n  —  m)  (8) 

where  f(m)  are  the  scaling  function  coefficients  of 
/(z)  at  scale  m,  4>{m)  are  the  corresponding  wavelet 
coefficients,  H(m)  and  G(m)  are  the  wavelet  trans¬ 
form  analysis  operators,  and  the  notation  *  indicates 
adjoint.  Finally,  Equations  7  and  8  indicates  that 
the  wavelet  coefficient  sequence  is  a  zero  mean  white, 
Gaussian  process  with  covariance  at  scale  m  given 
by  R(m).  This  particular  a  priori  model  defines  a 
stochastic  process  recursively  in  the  scale  parameter 
m  on  a  lattice  of  points  [1].  More  general  lattice-type 
processes  can  be  constructed  [1]  which  may  be  of  use 
in  future  work  on  multiscale  inversion  techniques. 

4.2  Inversion  Algorithms 

Given  this  multiscale,  stochastic  regularization  we 
may  immediately  construct  the  Maximum  a  posteriori 
estimator  for  the  wavelet  coefficients  of  f(z)  as  [6] 

4>map  =  A^©r[©A^©r  +  A^]_1j7  (9) 

with  A*  =  cov(x)  and  where  we  have  tacitly  assumed 
i  =  1.  The  form  of  the  gain  matrix  multiplying  rj  for 
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Figure  4:  Grayscale  Image  of  Absolute  Value  of  Gain 
Matrix  in  Equation  10  Using  Haar  Wavelet 


the  kernel  7\(z)  is  given  in  Figure  4  demonstrating 
that  the  stochastic  inverse  can  be  a  structured  and 
sparse  matrix  which  should  be  constructible  in  an  ef¬ 
ficient  manner. 

Construction  of  the  inverse  according  to  Equation 
9  is  not  efficient  due  to  the  necessity  of  directly  in¬ 
verting  the  matrix  in  brackets.  In  the  presentation 
accompanying  this  paper,  we  discuss  inversion  algo¬ 
rithms  which  utilize  the  system  defined  by  Equations  5 
through  8  to  synthesize  state-space-like  models  which 
allow  the  inversion  to  proceed  in  a  scale-recursive  man¬ 
ner.  The  model  so  constructed  must  meet  two  criteria. 
On  the  one  hand,  we  require  some  level  of  fidelity  from 
the  model  in  that  it  should  reproduce  the  action  of  the 
original  operator  to  within  acceptable  and  quantifiable 
limits.  On  the  other  hand,  the  model  must  be  of  suf¬ 
ficiently  low  complexity  so  that  it  lends  itself  well  to 
computationally  efficient  inversion  algorithms. 

Using  these  approximate  models,  we  identify  three 
broad  classes  of  inversion  procedures  which  may  be 
considered.  First,  the  models  may  be  both  suffi¬ 
ciently  accurate  and  of  sufficiently  low  complexity 
such  that  estimates  obtained  using  these  models  rep¬ 
resent  the  final  solution  to  the  inverse  problem.  Al¬ 
ternatively,  we  may  utilize  the  approximate  models  as 
pre-conditioners  for  other  solution  techniques  such  as 
successive  over  relaxation.  Finally,  it  may  be  possi¬ 
ble  to  improve  the  estimates  generated  under  out  ap¬ 
proximate  models  using  an  iterative  correction  scheme 


built  around  the  multiscale  model  itself  rather  than  an 
SOR-type  approach. 

5  CONCLUSION  and  FUTURE  WORK 

This  paper  has  presented  the  framework  for  a  mul¬ 
tiscale  approach  to  the  solution  of  one  dimensional  in¬ 
verse  problems.  We  have  discussed  the  utility  of  this 
approach,  demonstrated  its  potential  through  the  use 
of  a  relevant  example,  and  developed  the  basics  of  a 
multiscale  stochastic  regularization  method  for  solv¬ 
ing  the  problem.  In  the  presentation  accompanying 
this  paper,  we  discuss  scale-recursive,  computationally 
efficient  reconstruction  algorithms  derived  from  state- 
space-like  models  constructed  to  approximate  the  orig¬ 
inal  operator. 
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